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Topological phases supporting non-abelian anyonic excitations have been proposed as candidates 
for topological quantum computation. In this paper, we study disordered non-abelian anyonic chains 
based on the quantum groups SU{2)k, a hierarchy that includes the v — 5/2 FQH state and the 
proposed v — 12/5 Fibonacci state, among others. We find that for odd k these anyonic chains realize 
infinite randomness critical phases in the same universality class as the Sk permutation symmetric 
multi-critical points of Damle and Huse (Phys. Rev. Lett. 89, 277203 (2002)). Indeed, we show that 
the pertinent subspace of these anyonic chains actually sits inside the C Sk symmetric sector of 
the Damle-Huse model, and this Zfc symmetry stabilizes the phase. 



I. INTRODUCTION 

One of the major advances in the understanding of 
strongly correlated quantum systems has been the explo- 
ration of topological phases of matter. Originating with 
the discovery of the = 1/3 fractional quantum Hall ef- 
fect, topological phases have received much renewed in- 
terest with their recent proposed application to quantum 
computation'^. Under this proposal, quantum computa- 
tion is carried out by the braiding of the non-abelian 
quasi-particle excitations of the topological phase. The 
topologically protected degenerate space of ground states 
of the non-abelions serves as the memory, and the braid- 
ing induces unitary transformations within this Hilbert 
space^. The remarkable feature of this scheme is that 
the dimension of this space for N anyons grows as , 
with the quantum dimension d in general non-integer. 
This should be contrasted, for example, with N spin-1/2 
quasi-particles, whose Hilbert space has 2^ states. The 
non-integral nature of d reflects a unique non-locality 
of the anyon Hilbert space and makes decoherence-free 
quantum computation possible: no local perturbations 
can give rise to decoherence. 

The potential applications of topological phases with 
non-abelian anyons, and perhaps even more so their re- 
markable properties, warrant an analysis of interacting 
many-anyon systems. Possibly the most natural start- 
ing point is a 1-dimensional anyonic chain, analogous to 
regular spin chains such as the SU (2) Heisenberg model. 
Though such anyonic chains do not have truly local de- 
grees of freedom, they can be written as local systems 
with local constraints. It was found that even rather sim- 
ple to construct translationally invariant chains, based 
on SU{2)3 anyons, have an intricate structure, with the 
low energy degrees of freedom organized into the gap- 
less spectra of the c — 4/5 and c = 7/10 conformal field 
theories^. These chains interact via nearest neighbor cou- 
plings and are exactly solvable, though more generalized 
models exhibit an even richer set of critical and gapped 
phases'*. 



Whereas translationally invariant chains are the req- 
uisite first step, a more likely physical realization of a 
non-abelian chain, for example in a quantum Hall state, 
will have a strong degree of quenched disorder. Such 
chains are the focus of this paper. Even with garden- 
variety spins, the addition of disorder dramatically af- 
fects the physics, with the low energy behavior controlled 
by so-called infinite randomness fixed points^'^. A pro- 
totypical example is the random singlet phase describ- 
ing the state of a disordered spin-1/2 Heisenberg chain. 
Here spins pair up and form singlets in a random fashion, 
with most connecting near neighbors but some being very 
long ranged. This unique structure of the ground state 
leads to some rather unexpected or unusual properties, 
including algebraically decaying average correlations and 
energy-length scaling | In i?] ~ L'^ as opposed to the usual 
1/E ~ characteristic of pure systems. Furthermore, 
disordered spin-5' chains with 5' > 1/2 were shown to 
exhibit infinite-randomness critical fixed points with the 
universal exponent — 1/(25), and with the spin-state 
described as a competition between 2S + 1 domains, as 
we describe below. These fixed points were dubbed the 
Sn permutation symmetric points, with n — 2S + 1 the 
number of competing domains'"^ 
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Because of the unique structure of their Hilbert spaces, 
disordered anyonic chains are particularly amenable to 
treatment via strong randomness renormalization group 
methods, and indeed have been shown to exhibit infinite 
randomness fixed points*^'^^. They are thus an especially 
fertile ground for trying to discover and classify new uni- 
versality classes of strongly random behavior. Indeed, 
even though no new universality classes were found, Ref. 
[12] found a rich phase diagram for the SU{2)3, or Fi- 
bonacci anyonic chain, with a random singlet fixed point 
that can be destabilized by the addition of couplings fa- 
voring fusion into a non-trivial topological charge, and a 
resulting flow to a more intricate ^3 permutation sym- 
metric fixed point. It is notable that this 6*3 symmetric 
point is actually a critical phase, with no relevant per- 
turbations - this contrasts with the domain realization of 
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the Sn symmetric fixed point, which has n — 1 relevant 
perturbations^. 

The non-abelian models that are the subject of this 
paper are based on SU{2)k anyons with k odd. The level 
k signifies a truncation of all representations with spin 
S > k/2 - the anyons of SU{2)k simply correspond to 
the first k+1 irreducible representations of SU{2). This 
fact potentially suggests a subtle analogy between the 
SU{2)k systems and spins S truncated at k/2. In par- 
ticular, it raises the possibility of a relation between the 
Sk permutation symmetric fixed points of regular spin- 
(fc — l)/2 chains, and the infinite randomness phases aris- 
ing in SU{2)k anyons. 

In this paper we show that, indeed, for all odd fc, 
SU{2)k anyonic chains realize Sk symmetric infinite ran- 
domness critical phases. Crucial to our analysis is ex- 
pressing the 2-anyon interaction terms in a novel basis, 
one which behaves better than the standard projector ba- 
sis with respect to the renormalization group (RG) dec- 
imations. With this insight, we are able to solve the 
model, and in fact construct an explicit equivalence be- 
tween it and the C Sk symmetric sector of the domain 
model of Damle and Huse'', where the order 2k dihedral 
group Dfe in particular contains Z^. The symmetry is 
what stabilizes the phase, eliminating the k—1 relevant 
perturbations of the domain model. This phenomenon 
of a /c multi-critical point being stabilized by additional 
symmetry in the SU{2)k model is tantalizingly close to 
that discovered in [13] for the imiform SU{2)k chain. 
There, stability of an a priori k multi-critical point is 
guaranteed by an extra "topological" symmetry in the 
quantum system. It is conceivable that these two phe- 
nomena are closely related. 

The rest of this paper is structured as follows. In sec- 
tion II we briefly review necessary background on the 
strong randomness renormalization group procedure, as 
well as set up the construction of the anyon chain Hilbert 
space and Hamiltonian. In Section III we analyze the 
disordered anyon chain, introducing the novel basis for 
the interaction terms. We write down and solve the flow 
equations, finding a fixed point of the RG. In Section IV 
we construct an explicit equivalence between the disor- 
dered anyon chain and the ©fe symmetric sector of Damle 
Huse domain wall model, one that relates their respec- 
tive fixed points. We conclude in Section V, and relegate 
some technical derivations to the appendix. 



II. BACKGROUND AND SETUP 

A. Real Space Renormalization Group 

To find the ground state of disordered spin chains, 
Ma and Dasgupta introduced the strong disorder real- 
space renormalization group method^^'^^. The random 
spin- 1/2 Heisenberg model provides the simplest exam- 
ple amenable to such a treatment. The model is given 



by: 

■ff = ^ ■^i.i-l-l Si ■ Sj+i, (1) 
i 

where the couplings Ji^i+i > are positive and randomly 
distributed. Note that, as far as the Hilbert space is 
concerned, we have for two neighboring sites 

i®i=0©l, (2) 

and the interaction Sj • S^+i simply gives an energy split- 
ting between the singlet and triplet. The procedure now 
is to pick the largest Ji,i+i, which effectively decimates 
the excited triplet and leaves the ground state in a singlet, 
and do perturbation theory around that state. Quantum 
fluctuations then induce an effective coupling according 
to the so-called Ma-Dasgupta rule^^'^^: 

'Ji-l,i+2 = ^ 

So sites i and i + 1 are decimated and replaced with an 
effective interaction between i — 1 and i + 2. Iteration 
of this procedure produces singlet bonds on all length 
scales. This is the random singlet ground state. 

A quantitative description is obtained by tracking the 
RG flow of the coupling distribution. It is useful to em- 
ploy logarithmic couplings^: 

A,i+i=ln-^ (4) 

where ^ = maxjjj_j+i. In these variables the Ma- 
Dasgupta rule (3) reads 

A-l,i+2 = + Pi+l,i+2 (5) 

(up to an additive constant of In 2, which can be safely 
neglected). As the couplings get decimated Q, decreases. 
It is convenient to define the RG flow parameter as 

r = .n^, (6) 

where ilo is the maximal coupling of the bare Hamilto- 
nian. Let fr(/3) be the distribution of couplings. We can 
derive a flow equation for Priji) by decimating the cou- 
plings in the infinitesimal interval (3 = [0, dF] and seeing 
how their probabilistic weight is redistributed. We obtain 

P(0) dp, d(32S{P - /3i - (3)Pr{f3i)Prm 

The first term comes from the overall change of scale, and 
the second from the Ma-Dasgupta rule. These equations 
have a solution 

Prif}) = fe-f"/^ (8) 
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FIG. 1: a) Valence Bond Solid of type a = 2. b) Sample 
configuration at the infinite randomness fixed point. All the 
domains occur equally. The dashed lines represent domain 
walls, which contain the unpaired spins that form the low 
energy degrees of freedom. 

which is an attractive fixed point to essentially all phys- 
ical initial configurations. This solution permits us to 
read off features of the random singlet phase; for example 
one can with a little more work derive the energy-length 
scaling relation: 

Li/2_r = ln (!/£;). (9) 

which thus has the exponent: 

^ = 1/2. (10) 

The random singlet phase describes the universal low- 
energy behavior of several known one dimensional sys- 
tems, making it interesting to attempt to classify all such 
universal low energy fixed points of strongly random sys- 
tems in one dimension. Thus far all known universality 
classes are realized in the Damle-Huse hierarchy of per- 
mutation symmetric critical points^, indexed by a posi- 
tive integer n. In the construction of Damle and Huse, 
the system indexed by n is realized in a spin-n/2 SU{2) 
invariant Heisenberg model, with the sites represented as 
symmetrized tensor products of spin 1/2 's. The interac- 
tion Hamiltonian is 

H = Y,JiSi-Si+i (11) 

i 

where the couplings Jj contain dimerization 6 and ran- 
domness of strength R 

Ji = J[l + 5{-iy]exp{Rr]i). (12) 

Here the rji are random variables. Depending on R 
and 6 this Hamiltonian can realize a plethora of phases, 
which can be qualitatively understood in terms of Va- 
lence Bond Solids (VBS). In this picture, a VBS of type 
a € {0, 1, . . . , n} is constructed by pairing up a spin-l/2's 
into singlets over each even bond and n — a spin-l/2's 
over each odd bond (figure 1. a). This exhausts the spin 
degrees of freedom and defines a unique state. 



The Sn permutation symmetric infinite randomness 
fixed point is now realized multi-critical point in 

which all of these domains occur equally (figure 1. b). 
The domain walls contain unpaired spins: between a and 
a' we have a spin of magnitude \a — a'\/2. These spins 
interact via effective couplings whose magnitude is highly 
dependent on the domain that separates them, and whose 
sign is dictated by a consistency condition. Specifically, 
for three domains (Ti , (T2 , and C73 , the interaction between 
|fi — cr2|/2 and |a2 — 1J3I/2 is anti- ferromagnetic if ai —a^ 
and fT3 — fT2 are of the same sign and ferromagnetic oth- 
erwise. The fixed point turns out to contain an entirely 
random distribution of domains, described by a stochas- 
tic transfer matrix with all nonzero entries equal to 1 /n, 
an energy-length infinite randomness scaling exponent 
'ip = I In, and a logarithmic distribution of the coupHng 
strengths: P(/3) = (n/F) e""''/'^. 

We will see in the remainder of the paper that this 
multi-critical point will be realized (for n + 1 odd) as 
a stable phase of the SU{2)k anyon chain. First, how- 
ever, we need to review some background on anyonic spin 
chains. 



B. Hilbert Space and Hamiltonian of SU{2)k Anyon 
Chains 

A crucial part of our analysis relies on the specific 
properties of SU{2)k. The nontrivial anyons in this case 
correspond to the nontrivial representations of SU{2) 
at level k. There are k of these, labeled by their spin: 
0, 1/2, 1, . . . , k/2. The fusion rules are: 

min {i+j,k—i-j) 

j = ^ m, (13) 

m=\i-j\ 

where m is summing over the integers if i—j is an integer, 
and over the half integers otherwise. More information is 
encoded in the so-called F-matrix or set of q-G-j symbols 
of SU{2)k. Given 3 anyons ji, j2, and J3, we can either 
fuse ji and ^2 first into ji2 and then with into j, or 
we could first fuse ,72 and ,73 into ,723 and then with j\ 
into j. Both of these procedures generate a basis for the 
the Hilbert space of ground states of j\, j2, and J3. The 
transformation between these two bases is encoded in the 
i^-matrix iFf^JjU. The F-matrix of SU{2)k is written 
down in the appendix. 

Let us now construct the Hilbert space for the prob- 
lem. We have a chain {ip} of anyons indexed by an 
integer position p. From now on we will deal only 
with odd level k and integer values of the "spin" - thus 
ip G {1, . . . ,{k — l)/2}. We can do this because the inte- 
gers form a closed fusion sub-algebra of SU{2)k. Indeed, 
there is a Z2 symmetry that relates charge j and k/2 — j, 
and this symmetry will become important when we relate 
the anyon model to the Damle-Huse model. 

There are two equivalent ways to define the Hilbert 
space^^. The simplest construction is to label each "link" 
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FIG. 2: Sample link basis vector for an 5(7(2)5 chain. The 
anyons are represented by the black boxes at the bottom and 
labeled with their topological charge, 1 or 2. The links, or 
bonds, between the anyons are labeled with 0, 1, or 2, in such 
a way that the fusion rules at the trivalent vertices are obeyed. 



FIG. 4: Graphical representation of a projector onto total 
topological charge / acting on sites 2 and 3. The oval labeled 
by G represents a graph as described in the graphical defini- 
tion of the Hilbert space. The normalization constant C is 
chosen so as to make V a projector. 
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FIG. 3: F-matrix reconnection rules. The first is the so- 
called "no tadpole" rule, stating that any graph that can be 
disconnected by the removal of one edge is 0. The second 
states that removing a loop with label m is equivalent to 
multiplying by the quantum dimension d,n = (F™™)o. The 
third states that by performing the indicated reconnection 
on any local portion of the graph, we have the stated linear 
relation among the resulting graphs. 

between site p and p+ 1 with an integer anyon type lp,p+i , 
subject to the constraint that the fusion rules be obeyed 
at each site, i.e., 

^p,p+i ^ {l^p— i,p ip|,...,niin (lp—i^p-\-ip^k Ip—i^p ^p)} 

(14) 

See, for example, figure 2. The set of all such admissi- 
ble labclings defines a basis for the Hilbort space, which 
is simply the space of degenerate ground states of this 
configuration of anyons. 

We can impose the link basis constraints as high en- 
ergy 2-body interactions. The advantage of the link basis 
is then that it gives us a local way to describe the degrees 
of freedom in the problem. Indeed, we will see that the 
Hamiltonian defined below will consist of 3-local inter- 
actions. In order to define the Hamiltonian, however, it 
is useful to first describe a second, more abstract way to 
define the Hilbert space. 

Let us for convenience suppose the chain is finite, con- 
sisting of TV anyons. We first consider the space of all 
trivalent graphs, with endpoints on the N anyons, whose 
edges arc labeled by the non-trivial integer anyon types, 
and whose vertices satisfy the fusion rules. We take the 
Hilbert space generated by such graphs (modulo graph 
isomorphism) and quotient out the subspace generated 
by the f -matrix relations, interpreted as local reconnec- 



tion rules (see figure 3). To relate this graphical picture 
to the link basis, note that the link basis states can be 
viewed as labeled trivalent graphs, and that any other la- 
beled trivalent graph can be reduced to a superposition 
of these using F-matrix reconnection rules. The inner 
product of two graphs is defined by reflecting one of the 
graphs and concatenating it with the other along the N 
nodes. 

An advantage of this graphical picture of the Hilbert 
space is that it makes it easy to describe the interaction 
terms occurring in the Hamiltonian. The Hamiltonian 
is a sum over p of pairwise interactions between sites p 
and p + 1. These are constrained by the SU{2)k sym- 
metry, and therefore a linear combination of projection 
operators onto some total topological charge /: 

n = Y.JL+i'Pp2+i (15) 

These projection operators V have a graphical represen- 
tation, and in the abstract graph basis their action on a 
particular graph G is simply given by concatenation of 
V with G, up to a normalization constant (see figure 4). 
Their action in the link basis can be worked out by con- 
catenating V with a particular link graph and then using 
F-matrix rules to reduce the resulting graph to a linear 
combination of link graphs. In this manner it is apparent 
that the action of Vp^p+i depends only on Ip-i^p, lp,p+i, 
and Zp+i,p+2; it is thus a 3-local operator. 

III. ANALYSIS OF THE DISORDERED SU{2)k 
CHAIN 

A. Convenient Basis for the RG 

We would like to apply the real space RG procedure 
to the disordered SU{2)k chain (15) in hopes of find- 
ing infinite randomness fixed points. Let us first review 
what happened in our previous analysis [12] of the case 
fc = 3, i.e., the Fibonacci chain. S'C/(2)3 contains one 
nontrivial anyon of integer charge, the so-called t anyon, 
and the chain is simply an array of these. In the strong 
randomness limit, we applied the Ma-Dasgupta rule to 
the strongest bond, which was a projection operator on 
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FIG. 5: Graphical definition of tfic operator Q^p p^i- 
tliat it is nonzero only when the fusion rules at the (ip,; 
and vertices are obeyed. 
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a pair of neighboring r's. These r's could fuse to one 
of two possible states, either one with trivial total topo- 
logical charge or another r {t ® t = © r), and the 
bond projected onto one of these, leading to either the 
elimination of both anyons, or their merger into one. In 
either case, we were left with effectively another realiza- 
tion of the Fibonacci chain, with either 1 or 2 fewer sites, 
allowing us to iterate the procedure. 

For SU {2)h a new complication arises. This time, when 
we pick the largest bond to decimate, the generic situa- 
tion is that there arc more than two possible fusion prod- 
ucts for the corresponding pair of anyons. For example, 
in SU{2)5, the fusion rules are: 



101 = 0©!© 2, 

202 = 0© 1, 
20 1 = 1 ©2. 



(16) 



The first rule, regarding the spin 1 representation, con- 
tains 3 possible fusion products on the right hand side. 
We would like to be able to keep just the lowest energy 
of these, in order to merge the two anyons into a sin- 
gle new effective anyonic site. In general, however, we 
are not allowed to do this: when there arc more than 
two fusion products, there will be more than one energy 
splitting associated with the bond. While we can dec- 
imate away the largest one, there may be coiiplings on 
other bonds that need to be decimated before the smaller 
couplings on the original strong bond. This decimation 
of only one fusion product leads to a situation where on 
the one hand we need to enforce a constraint on the two 
anyons, but on the other we are not allowed to merge 
them into a single effective anyon. This impasse makes it 
seemingly impossible to carry out an iterative real space 
RG analysis. 1^ 

In the rest of this paper we demonstrate the existence 
and stability of an infinite randomness fixed point for 
SU{2)k chains that circumvents the above difficulty. The 
idea is to construct a Hamiltonian out of two-site opera- 
tors Qp,p+i 



p+1 Qp,p+i 



(17) 



such that H retains its form under an RG procedure 

where we truncate all excited fusion products. In other 
words, the effective operators generated from first and 



second order decimations are all proportional to Q. The 
a priori assumption that all such excited states can be 
truncated, which is not valid in general, is then justified 
in this particular case if we can show that the resulting 
RG leads to a strongly disordered set of couplings Jp,p+i- 
This is because when we express the operator Qp,p+i as 
a linear combination of projection operators 



(18) 



the differences between the coefficients c / remain of order 
1, independently of the broadness of the distribution of 
log Jp_p-|_i. Hence the energy splittings for each bond are 
all of the same order, and with very high probability all 
get decimated in one fell swoop in the large disorder limit. 
The assumption that all excited states can be truncated 
is then justified and the scheme is self-consistent. Of 
course, this does not rule out the possiblity of more exotic 
phases where the interactions are not built out of only the 
Q operators - these phases, if they existed, would not be 
amenable to treatment by this method. 

We now claim that the correct operators to tisc arc the 
ones defined graphically in figure 5 and denoted Q^^\ 
They can be thought of as an exchange of an anyon of 
topological charge /; their action on any graph G is sim- 
ply by concatenation with G, as for any graphically de- 
fined operator. The Q^^^ differ from the projection op- 



erators V. 



(9) 

P,P+1 



by an .F-matrix move, and can of course 



be expanded as linear combinations of the Vp p_^_i. They 
all have well defined scaling under the RG, and, as we 
show in the appendix and discuss in more detail below, 
the most relevant one is Q(i). Thus from now on we will 
consider a Hamiltonian of the form 



H = y^^Jp,p+iQ. 



(1) 

■p,p-i-i 



(19) 



where the couplings Jp,p+i are disordered. In the ap- 
pendix, we show 



min (i-\-j^k—i—j) 
f=\i-3\ 



(20) 



where 



m = [fii + [f+^il-i\^'-j\il 



+ (21) 



is an increasing function of /. Here i and j arc the topo- 
logical charges of the two neighboring anyons on which 
Q^^^ acts, and the g-numbers are defined as 



(22) 



with q = e'^'/('=+2). 

Thus, depending on the sign of the coupling Jp.p+i, 
the anyons at p and p+1 can fuse to either an anyon 
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FIG. 6: Graphical representation of the effective operator be- 
tween the composite i—j and its neighbor m, generated at first 
order in perturbation theory. It is equal to ■p}2~^^Q23'^i2 
By performing F-matrix manipulations within the dashed 
box, we see that this effective operator is just a multiple of 

■^comp. 3* 

of charge \ip — ip+i\ or one of charge min [ip + ip+i, k — 
ip — ip+i)- As discussed above, we decimate all the other 
fusion products. This scheme will be self consistent if 
wc show that the system flows to strong disorder, i.e., 
that the distribution of the log Jp,p+i broadens out. To 
analyze the flow, we first need to work out the decimation 
rules. 

First, let us focus on first order decimations. Here 
the strongest bond Jp,p+i fuses the p and p + 1 anyons 
into a composite, and the effective interactions between 
the composite and its neighbors depend on the original 
< p—l,p > and < p+l,p+2 > interactions. Now, for the 
Hamiltonian (19) we see from figure 6 that these effec- 
tive interactions are proportional to Qp-i^p and Qp+i.p+2 
- the figure represents graphically the first order pertur- 
bation calculation^^. In fact, wc can say more: figure 6 
makes it clear that not only are the Q*^^) preserved, but 
so are all the Q^^\ More precisely, if we think of a first 
order decimation as a linear mapping from the space of 
interactions between say p+1 and p+2 {or p and p—1) 
to the space of interactions between the composite and 
p + 2 (or p — 1), figure 6 makes it clear that the operators 
ai'e eigenvectors of this mapping. 

Second order decimations follow a similar paradigm, 
and in the appendix we show that, again, the Q^^^ arc 
eigenvectors of the second order decimations, and derive 
the real-space Ma-Dasgupta decimation step, which (for 
Q(i)) reads: 

neff = Q ^ Q[% (23) 

^2,3 

as in Eq. (51), with Cj given by the expression (52). 
For both first and second order decimations, we also cal- 
culate in the appendix the corresponding eigenvalues for 

all P, and show that in each case they are maximized for 
(3=1. Thus Q^^-* is the most relevant operator, which 
is why we chose to construct the Hamiltonian out of it 
in (19). ft is stable with respect to perturbations by the 

for (3 > 2, which, having smaller eigenvalues un- 



der decimation, arc irrelevant. 

In fact, the stability argument is not entirely rigorous, 
but follows the usual line of justification for the validity of 

strong randomness RG. Basically one can show that the 
addition of a small amount of Q^'^p+i for /3 > 1 doesn't 

change the decimation rules, and so these Qp'^^+i decay 
away under the RG, up to bad spots, or "cancers", that 
occur with frequency that vanishes with increasing ran- 
domness. We can invoke the standard line of reasoning 
used to justify strong randomness RG in the first place 
to argue that these do not destabilize the fixed point, 
though ultimately this should be decided by numerical 
simulation. 

Thus, with the ansatz (19) for the Hamiltonian, we 
have a consistent framework for the RG that elimi- 
nates the potential multitude of widely distributed en- 
ergy scales associated to each bond. Instead, we have 
only one energy scale for each bond, multiplying the oper- 
ator Q^ppj^i- In the strong randomness limit, decimation 
of the strongest bond results, with probability approach- 
ing 1 , in the decimation of all its excited states, leaving 
a or 1 anyons in place of 2. This decimation process 
preserves the form of the interactions Q^pp^i- All that we 
have left to do is to show that under the RG the Hamil- 
tonian (19) does indeed flow to strong randomess. This 
we do in the next subsection. 



B. Flow Equations 

Let us now solve the model (19). To do this, wc need 
to derive the flow equations that describe real space RG. 
We are dealing with an ensemble of chains, where not 
only the coupling strengths and signs, but also the anyon 
types arc chosen according to some probability distribu- 
tion (sec also [8]). We make the ansatz that the cou- 
pling strengths, signs, and anyon types are completely 
uncorrelated from each other, and uncorrelated among 
the different sites/bonds. Also, for simplicity wc ana- 
lyze only the case in which the bond strength probability 
distribution is symmetric with respect to sign change. 
The stability analysis in this framework might in prin- 
ciple miss asymmetric relevant perturbations, and in- 
deed other non-independent distribution perturbations, 
but the exact mapping to (a symmetric subspace of) 
the Damle-Huse model discussed in the next subsection 
will show that there are none, at least at the level of 
analysis in [7]. 

Let R{i), i = 1, . . . , (A: — l)/2 be the probability distri- 
bution for the (integer) anyon types, and P(/3) the log- 
arithmic bond-strength probability distribution, normal- 
ized to integrate out to 1 /2 because of the two possibili- 
ties for the sign of the coupling. Here f3 = log(Cl/ J) is the 
logarithmic coupling, and fl the energy cutoff. Consid- 
ering both first and second order decimations, we obtain 
the following infinitesimal transformation for the joint 
probability distribution R{i) P{P): 
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R{i) P{(3) Rit) P{(3 + rfr) + drP(O) {R ® R){i}P{P) + rfrP(O) {R^R){Q) R{i) {P (g) (24) 

The notation is defined below. Here the first term comes from the cutoff rescaling, the second from first order 
decimations, and the third from second order decimations. Equation (24) is equivalent to the following two norm 
preserving transformations of R{i) and P{/3)'- 

R{i) Rii) + P(0) dr {R «) R){i) - P(0) dT (1 - (i? o i?)(0)) i?(i) , . 

P{(3) P(/3) + (TP' {(5) + P(0) dV {R (8> P)(0) (P (8> P)(/3) - P(0) dV (1 - (P ® P)(0)) P(/3) ^ ^ 



The notation is as follows. We define the convolution 

(P ® P) (/3) = 2 dfi' P{(3') P{(3 - (3') (26) 
Jo 

with an extra factor of 2 to account for normalization of 
P, and we let 

(k-l)/2 

{R^R){i)= v{i,j,l)R{j)R{l) (27) 
I 



A solution to these equations is 

Let us analyze the stability of the solution in Eq. (29) . 
First let us look at R{i). Consider a perturbation of the 
form R{i) = Ro{i) + €i. Using the easily derived fact that 
to linear order {R (g) R){i) = (Po ^ Ro)ii) — 2ej and the 
fact that J2i = we get the RG flow of Cj: 



^ = -P(0)e, (29) 

which, since P(0) is always positive, shows that Ci always 
decays. Now consider P(/3). At linear order the variation 
in (P P) (0) vanishes, so the analysis of the stability of 
P(/3) is as in all the other examples of strong randomness 
RG where this solution occurs^'^. As mentioned before, 
wc don't consider perturbations asymmetric with respect 
to the sign of the coupling, but we expect these to be 
stable as well by arguments similar to those in [12]. We 
also show now via explicit mapping to the Damle-Huse 
model that they indeed are stable. 



I 

where ry(i, j, is equal to 1 if i = |j — Z| or i = min {j + 

l,k — j — I) and is otherwise (the two possibilities cor- 
respond to the two possible signs of the coupling). 



We can also write down the integro-differential flow 
equations corresponding to these infinitesimal transfor- 
mations (the r dependence is implicit): 



(28) 

I 

IV. RELATION TO THE DAMLE-HUSE 
DOMAIN WALL MODEL 

As mentioned above, the fixed points wc found above 
must somehow be related to the Damle-Huse fixed points 
of abelian spin chains. In this section we present and 
discuss the mapping between the Damlc-Husc domain 
model with n = k domains and spin S = {k — l)/2, and 
our SU{2)k anyonic chains, and show that this mapping 
gives an equivalence between the permutation symmetric 
multi-critical point of the domain model and the fixed 
point of the non-abelian anyonic chain. In addition, we 
show that the anyonic fixed point is actually a stable 
phase. 



A. Mapping Between the Two Models 

Some intuition for the mapping between the spin-S* 
Damle-Huse chain with S = (fc — l)/2 and the SU{2)k 
anyon chain can be obtained by inspecting the Brat- 
teli diagram in Fig. 8. Naively, we like to think of 
the possible topological charges i = 0, 1/2, fc/2 of 
the SU{2)k tensor category as somehow related to spin- 
i representations of SU{2). This naive notion is not 



£^ = P(0) (P ® P)(i) - P(0) (1 - (P P)(0)) R{i) 
^ = P'iP) + 2P(0) (P ® P)(0) dP' P{I3') P(6 - 6') - P(0) (1 - (P f: P)(0)) P{6) 
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a=3 a=i a=4 a=2 a=i a=3 a=i a=2 a=i 

I I domain model i i i 

— decimation ^ m- ^ — ^ 




I +1 I +1 I -1 I 

s=1 s=3/2 s=1 s=1/2 



1+1 I +1 I 
s=1 s=1/2 s=1/2 




anyon model 
decimation 



+1 



-1 +1 



FIG. 7: An example illustrating commutativity of the equivalence mapping between the spin-2 Damle-Huse chain and the 
SU{2)^ anyon chain and (first order) decimations. A spin-2 configuration has five different domains, while the anyon model has 
two non-trivial charges: 1 and 2. The domain model configuration at the upper left, has domains labeled by o" = 3, 1, 4, 2, 1, signs 
of couplings -1-1,-1-1, —1, and strongest bond spanning the 4 domain. These are mapped to the SU{2)z domain (down arrow), 
with topological charges following (31) as j( """il -j^ g^jjjj signs of the couplings following (32), s'(m) = (— l)('^»"+i~°"»"-i) . 

Decimating bond 4 (right arrow) in the anyon chain follows the principle stated in the appendix: when two anyons i and j fuse 
ferromagnetically to i-f j < fc/2, the effective couplings are the same as the original couplings (as opposed to other cases where 
the sign may flip). Starting again from the top left corner, carrying out a real space decimation (right arrow) in the Damle-Huse 
chain first, and then the mapping to the anyon model (down arrow) results in the same configuration, as illustrated. That 
this is true in general requires one to check a few more similar cases. We thus see that the mapping from the spin-(fc — l)/2 
Damle-Huse chain to the SU{2)k chain and real-space RG steps commute. 




FIG. 8: Bratelli diagram for 5(7(2)5. The mapping / defined 
below simply reflects the topological charges along the vertical 
axis for non-integer arguments, mapping non-integer charges 
to integers. Bars denote the cases where the argument was 
non-integer and hence the action of / non-trivial. 



quite correct, because of the special constraints that 
the Hilbert space truncation presents. As it turns out, 
both i and fc/2 — j essentially represent the same non- 
trivial topological charge. Therefore the distinct non- 
trivial topological charges can be indexed by integer z's: 
I = 1, 2, . . . , fc — 1/2; the half-integer values of i can be 
turned into integers through i — > fc/2 — i. This is, for ex- 
ample, the reason for our ability to restrict our rendition 
of the SU{2)5 fusion algebra to the rules in (13) using 
only the charges 1 and 2, alongside the trivial (vacuum) 
charge 0. 

The association of i = 0, 1/2, fc/2 with spin-i 
SU{2) representations provides the correct intuition for 
the mapping between the anyon model and the spin 
(fc— 1)/2 Damle Huse domain model. In the spin-(fc— 1)/2 
Damle-Huse chain, each site appears as a domain wall be- 
tween two domains, say am and cTm+i - in order to main- 
tain transparent notation we label such a site by the pair 



(to, m+1). The spin of site (m, m + 1) is expected to be 



S, 



m,m+l 



(30) 



The range of possible spins is 1/2 < S,n,m+i < {k — l)/2, 
just like i above, excluding the two trivial (vacuum) 
charges i ~ 0, fc/2. Thus in the mapping to the Damle 
Huse model, the natural thing to do is to identify the 
domain wall spins with the anyon charges. However, 
we would also like to restrict to the integer topological 
charges. This is naturally done through the mapping: 



fiS) = 



s 



S eZ 



s s- 



e Z 



(31) 



with 25 e Z, an integer, and < S < fc/2. So for 
example, for fc = 5, we have /(2) — /(1/2) = 2, and 
/(l) = /(3/2) = l. 

Now let us define the mapping a little more for- 
mally. A configuration in the Damle-Huse model is 
completely specified by a sequence {(Tm, Pm} of domains 
am GO,l,...,fc — 1 and log couplings f3m between them, 
since the signs s(m) = ±1 of the couplings are uniquely 
determined by this data. As described above, we map 
this configuration to the sequence {i(m,m+i} = /(Icm+i ~ 
tTm|/2)} of SU{2)k anyons, with the log of the couphng 



Jm between anyons i(^m-i,m) 
and the sign of Jm, given by 



and i{m,m+i) given by /?„ 



s'(m) = (-l)('^™+i 



^s(to). 



(32) 



This choice of sign will reproduce the prefered fusion 
channels upon mapping a spin-(fc — l)/2 chain (in its 
domain wall representation) to an anyonic problem. 
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This mapping commutes, by construction, with both 
the first and second order real-space RG decimations. 
Since this claim is the key point in the proof, for clar- 
ity, we illustrate it with a specific example in figure 7. 
In the caption we explain how we get the same anyon 
configuration and couplings, including signs, irrespective 
of whether we do the decimation in the domain model 
and then map to the anyon model, or first map to the 
anyon model and do the decimation there. This means 
that the two operations commute. Second order decima- 
tions (not illustrated) are even easier to handle. Here, in 
both the domain model and the anyon model, the sign of 
the coupling is given by the Ma-Dasgupta rule (3), and 
hence commute. We have thus defined a mapping from 
the configuration space of the Damle-Huse model to that 
of our anyon model, and this mapping respects the RG 
evolution. It is quite remarkable that despite the two 
different origins of the interaction couplings' signs (one 
from the Damle-Huse domain model rules, and one from 
an _F-matrix calculation) , the signs conspire to make the 
two operations commute. This is presumably a deep re- 
flection of the fact that the SU{2)k tensor categories were 
constructed from the spin-representation of 5(7(2). Also, 
note that while the Ma-Dasgupta rules in the two models 
may have different multiplicative constants, leading to a 
small difference in the logarithmic couplings between the 
two models, this difference is unimportant in the large 
disorder limit. 



B. Symmetry Considerations and Elimination of 
Relevant Perturbations 

Let us examine the properties of this mapping of con- 
figuration spaces. First of all, we claim the map is 
onto, i.e., given any configuration {i(m-i,m)j Pm, s' (m)}, 
there's a domain configuration {<7rm Pm, s(m)} that maps 
onto it. To see this, we first pick any < cri < fc— 1. Then 
we must pick < (T2 < k — 1 such that /(|(J2 — o"i |/2) = 
i(i.2)- It is easy to see that there are precisely 2 choices of 
such CT2 (naively we may think there are 4, given the 2 to 

1 nature of both / and the absolute value mapping, but 

2 of those choices arc not between and fc — 1). Now we 
must similarly choose 0-3; however, in this case we also 
have to choose it in such a way that the sign s(2) comes 
out correctly. This constraint uniquely determines (73, 
and in fact all the other (7,^ (for m > 2 and m < 0) are 
uniquely determined in this manner. 

Thus, we have not only shown that the mapping is 
onto, but also that each anyon model configuration has 
precisely 2k domain configurations that map to it, for 
the k choices of cti above, and the two choices of a2- 
Indeed, this 2fc-fold degeneracy is easy to understand. 
First of all, since the anyon charges are functions of only 
the differences between the am, we can add a constant 
c to all the am (mod) k without changing the anyon 
charges. In fact, this is somewhat subtle. For example, 
if we have am < o'm+i, the corresponding anyon charge 



is /(((Tm+i — am)/2). Now, if we add a constant c such 
that < am + c < k but am+i + c> k, i.e., so that am+i 
cycles through, the new anyon type is 

f{[{(Tm + C) - {am+1 + C - k)]/2) = 

fik/2 - {am+1 - Om)l2) = (33) 

/((Cm-l-l — 0'm)/2) 

by the definition of /. Thus the anyon type still remains 
invariant. The signs of the domain model couplings also 
get modified under such a cyclic shift of am+i, but this 
is precisely canceled by the corresponding modification 
of the s'{m) = (-l)('^"'-i-'^"'+i)s(m) sign rule. This cy- 
cling accounts for a fc-fold degeneracy; the factor of 2 
comes from the flip am —>■ k — am, which preserves ab- 
solute values of differences and signs of the couplings as 
well. Thus the set of 2k pre-images of any anyon model 
configuration is simply an orbit of the dihedral group Dfc , 
viewed as a subgroup of the symmetric group Sk acting 
on the configuration space of the Damle-Huse model. 

Now that we understand the nature of the mapping 
in Eq. (31) and (32), we can explicitly verify that the 
inverse image of the fixed point ensemble (29) under this 
mapping is precisely the Damle-Huse fixed point. This 
confirms that we have a map that identifies the two flxed 
points and commutes with RG evolution. Furthermore, it 
means that the physical properties of the two systems are 
the same, except those that might be affected by the 2k 
to 1 nature of the mapping. One of these is the existence 
of relevant perturbations - we will prove that the anyon 
model has no relevant perturbations, making it a stable 
phase. 

Before giving the formal argument regarding lack of 
relevant perturbations, we illustrate what happens with 
a physically appealing picture. Specifically, the domain 
model has fc — 1 relevant perturbations^ each of which 
can be described in terms of one domain falling out of 
favor with respect to the rest (there is a linear constraint 
since they can't all fall out of favor simultaneously). 
One might try to construct relevant perturbations of the 
anyon model by mapping these relevant perturbations of 
the domain model, as follows: each relevant perturbation 
can be thought of as a functional on the conflguration 
space of the domain model, so one can, for a given anyon 
model configuration, sum up the values of the relevant 
perturbation on all 2k of its pre-images (which are do- 
main model configurations). This sum, however, turns 
out to be 0, so no relevant perturbation in the anyon 
model can be constructed this way. 

Let us proceed by providing a formal proof based on 
symmetry. The proof is by contradiction. Supposing we 
had a relevant perturbation of the anyon model, we could 
then pull it back to a relevant perturbation of the domain 
model (given a map M : X —>■ Y oi spaces, the pull- 
back map on function spaces F{Y) F{X) is defined 
by / ^ 5 where g{x) = f{M(x))). Because of the 2k to 1 
nature of the mapping, this would yield a Dfc symmetric 
relevant perturbation of the domain model; in particular 
it would also be C Dfc symmetric. The fc — 1 relevant 
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perturbations of the domain model, however, are not 
symmetric, because they are described in terms of one 
of the k domains falhng out of favor with respect to the 
others. Indeed, they form a fc — 1 dimensional non-trivial 
irreducible representation of Zfc. Our putative relevant 
deformation is symmetric, i.e., lies in the trivial rep- 
resentation of Zfc. Thus we have found a non-existent 
relevant perturbation of the domain model, a contradic- 
tion. Of course, this analysis does not include possible 
perturbations by the addition of interactions Q^-^' with 
/ > 1, but we have already shown that the Hamiltonian 
is stable with respect to such perturbations in section III, 
i.e., we showed they arc irrelevant. We have therefore re- 
alized all of the odd k St symmetric multi-critical points 
as stable phases in the SU{2)k anyon chains. 



V. CONCLUSIONS 

In this paper we analyzed in full generality spin-chains 

made of non-abelian quasiparticlcs arising in the tensor 
categories of SU{2)k for odd k, in the limit of strong ran- 
domness. We have found a realization of the odd k Sk 
symmetric infinite randomness multi-critical fixed points 
of Damlc and Husc^ as critical stable phases of disordered 
SU{2)k anyon chains. We have shown that the SU{2)k 
fixed point is stable by analyzing the RG flow equations 
around it, and also by explicit mapping to the Damle- 
Huse model. Key in our analysis was our use of a basis 
of interaction operators Q*^/) that behave well with re- 
spect to real-space decimations, i.e., operators Q^/) that 
have well defined scaling under the RG. We found that 
Q*^l) is the most relevant operator, and thus the only one 
appearing in the effective Hamiltonian (19) at the fixed 
point. This effective elimination of all but one interac- 
tion operator is what resolves the a priori problem of 
having a potential multitude of energy sc;ales associated 
with the multiple fusion products of each neighboring 
pair of anyons. Indeed, at the fixed point (19) each bond 
is characterized by only one energy scale: the coefficient 
in front of the Q*^l) operator. 

Recall that the motivation for our study was the search 
for new universality classes of infinite randomness fixed 
points. In that sense, our analysis led to a disappoint- 
ment: The SU{2)k anyonic chains exhibited the same 
behavior as random spin-(fc — l)/2 chains, as though the 
differences between the extremely distinct Hilbert spaces 
of the two systems were essentially washed out in the 
strong randomness limit, leading to the same infinite- 
randomness fixed points. Nevertheless, a crucial differ- 
ence arose: the permutation symmetric fixed points mark 
stable phases of the SU{2)k random spin chains, as op- 
posed to unstable points in the ordinary spin chains. 

One natural question is then whether SU{N)k for 
higher A'' behaves any differently. It is plausible that 
using an approach similar to the one in this paper, with 
the relevant interactions as exchanges of a certain anyon 
type, will yield the already known fixed points, essen- 



tially because the charge of the exchanged anyon will pick 
out a preferred SU (2) C SU (N) and will decompose the 
problem to the SU{2) cases already studied. This does 
not, of course, rule out the possibility of more symmetric 
and exotic fixed points, which may arise with some fine 
tuning, for instance. 

Another line of investigation deals with the relation to 
the uniform SU (2)^ anyonic chains of 13. There, a ''topo- 
logical" symmetry stabilizes an otherwise k multi-critical 
point CFT low energy spectrum. A physical interpreta- 
tion of this phenomenon is given in terms of separating 
out the left and right moving modes of the chain, while 
creating a different topological liquid between them - the 
topological symmetry then eliminates relevant tunneling 
operators between the two modes. One could then add 
some very weak disorder to this system: at short dis- 
tances the picture of two modes separated by a topolog- 
ical liquid is preserved, while at long distances the disor- 
der grows and the dynamics is controlled by the infinite 
randomness fixed point discussed in this paper. It is in- 
teresting to try to find a physical picture for the infinite 
randomness phase - perhaps with the intervening liquid 
having broken up into disconnected islands - which may 
also yield the stability argument of the fixed points we 
found. 

On the other hand, a picture that is similar to the one 
of separating the right and left moving modes away from 
each other but which also applies to the strongly disor- 
dered system may provide clues to the understanding of 
the behavior of non-abelian anyons interacting on ran- 
dom planar graphs. Indeed, it would be interesting to 
see if any of the ideas developed in this paper have appli- 
cation to, say, a two (or higher) dimensional disordered 
lattice of anyons. So far we haven't made progress in this 
direction. 
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VI. APPENDIX 

In this appendix, we derive equation (20) for the energy 
spacing between the different fusion channels in the Q^^^ 
interaction, show that Q^^^ is the most relevant operator 
for both first and second order decimations, and derive 
the sign rules for the changes of the sign of the couplings 
under first order decimations. For the analysis we will 
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need the following expression^" for the F-matrix (or 6-j symbols) of SU{2)k: 



A/L2ii2 + MqV^hi + 1J9 A(ji , i2, ii2) A(ji2 , is, i) A(j2, is, i23) A(ii , i2s, i) 



(34) 



■j2 - jl2 J ,! Lz-J"i2 -33 -j] q ! Vz-j2 -h~h3\ q ' Vz-jl -J23 -jj g! 

X 1 

\jl +j2+j3+j-2j q ! LJI+J12+J3+J23 -«J (,! Lj2+jl2+j+j23 -^J g! 

I 



} 



Here the q-numbers are defined as in (22). The sum is 
over all z for which all the (jr-factorials arc well-defined, 
i.e., such that the arguments are > 0, and 



A(a,6,c) 



[-a + h + c\q\\_a-h + c\ 



[a + 6 + c + IJ 



9- 



Let us first analyze the operator Q 



(1) 



(35) 

To simplify 



notation, we denote {ip 



l} by {i, i} and drop the sub- 



scripts on Q^^^ and the projectors V'^^'' . Applying an 
i^-matrix move and noting the normalization on the pro- 



jector Vppj^i, we obtain 



Q(l) 



min {i-\-j,k—i 

E 

S=\i-3\ 



-j) 



(36) 



We would like to know the /-dependence of the coef- 
ficient in front of V^-^\ Plugging into (34), up to an 
/-independent prefactor, the coefScient is 



(L/J^ + L/ + iJ^-LN-j|J^-L* + j + iJ^) (37) 

giving us (20) as desired. Note that the quantity in the 
brackets is an increasing function of /. 

Let us now work out the decimation rules for the 
Hamiltonian (19), and in particular show that the Q 
operators are eigenvectors with respect to decimations, 
with Q(^) having the largest eigenvalue. Pick the largest 
coupling Jp,p+i. According to (37), anyons p and p + 1 
will be fused into an anyon of total topological charge 
\ip — ip+i\ for Jp,j,+i > and into an anyon of charge 

I 



max {ip + ip^i,k — ip — ip+i) for Jp,p+i < 0. Suppose 



first that this composite anyon charge is nonzero. In this 
case we need to use first order perturbation theory to 
work out the effective coupling of this composite anyon 
to its neighbors. For simplicity we only consider ip+2, the 
neighbor to the right, and to simplify notation we denote 
ip,ip+i,ip+2 by i,j,m. 

Let's first deal with a specific case, say having i > j 
fuse to i—j. According to first order perturbation theory, 
the effective coupling between i—j and m is given by 



(38) 



Graphically we can see (figure 6) that this effective cou- 
pling is still a multiple of Q^^\ There is a finite factor 
in front of the coefficient, as well as a possible sign, but 
we already see that the form (19) is preserved by first 
order decimations. In fact, it will be useful to go a little 
further and examine the effect of such a decimation on 



the other interactions Q, 



p-|-l,p-|-2- 



Again, we graphically 



see that the effective coupling will be a constant multiple 



of Q 



p+l,p+2- 



We are only interested in the sign and (3 



dependence of this constant, which turns out to be equal 
to the constant C in figure 6, this being the evaluation 
of the graph in the dashed box. This is 



(39) 



where a is the fusion product of i and j, in this case a = 
i — j. Plugging into (34), we see that the /3 dependence 
of (39) is 



M(/3) = (-1)^ (L2i - PU\[2i -2j + p + l},\[2j + P + l\,\[2i - 2j - /?],!)"'/' (40) 
We observe that the magnitude of this quantity decreases as a function of (3. To see this, note that 
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/x(/J)/m(/3 - 1) 



L2j^/3+lJq 

L2j + /?+lJ, 



1/2 



L2z-2.7-/?+lJ, 
L2i-2j + /?+lJ, 

I 



1/2 



(41) 



J Q 

2 2 



3^ 3 



J 



J 



FIG. 9: The setting for second order decimations. The anti- 
ferromagnetic coupUng J2 is much larger than Ji and J3, caus- 
ing the two middle anyons to combine into a charge singlet. 



It's easy to see using the explicit expression for the 5- 

numbcrs in terms of roots of unity that h{j) = h{k/2 — 
and that h{j) is an increasing function of j. There- 
fore 



h{j)h{i-j)<h{j)h{k/2-j) = l 



(42) 



Thus we have shown that Q^^' has the highest eigenvalue 
under first order decimations. We've only considered the 
case i > j, a = i — j, but similar arguments apply to 
other cases (for the case of a = A; — i — j we need to 
use the symmetries of the F-matrix discussed in [17]) to 
show that /u(/3) decreases as a function of (3. 

The one thing we need to know explicitly for the map- 
ping to the Damle-Huse model is whether first order dec- 
imations flip the sign of the neighboring couplings, when 
those are of the form Q^^^. For the case just consid- 
ered, i > j,a = i — j, the sign is (— I)'' = —1. For the 
rest of the cases the sign can be read off from the fac- 
tors of (—1)^ in (34): when i < j, a ^ j — i, the sign 
is -hi; when a = i + j < k/2, the sign is -|-1; and when 
a = k — i — j< k/2 the sign is —1. 

Consider now the second order decimations. The setup 
here is that we have four consecutive anyons, with the 
middle two anti-ferromagnetically fusing to the trivial 
channel, so we can label their topological charges a, j, j, b. 
The picture is as in figure 9. The bare Hamiltonian is 



Ho = J2{Qi'^-<Qi'^ 



>) 



where < 



> is the expectation value of 22^-* in the 



state where the fusion product of the two j's is trivial. 
The interaction Hamiltonian is 



H' = J^Qf + J3Ql 



(44) 



FIG. 10: Graphical expression for the effective interaction be- 
tween a and h generated at second order in perturbation the- 
ory. Here Ho is as defined in tiie text. Its inverse is computed 
and expressed in graphical form by expanding in projection 
operators. This is one of the two terms that contribute to an 
energy splitting between the fusion products of a and b. Note 
that both the original interaction with a and b exchange the 
same topological charge /3 - otherwise the second order cor- 
rection is zero by topological charge conservation, as is clear 
from the figure. 



H,s = rH' Hg^H' (45) 

Hes acts on the subspace where the two middle anyons, 
both of charge j, fuse to the trivial channel. The in- 
verse Hq^ is well defined in this expression because H' 
acting on this subspace gives a vector orthogonal to the 
subspace. We have, up to a multiple of the identity op- 
erator. 



H, 



J1J3 



eff 



(1) 



+ ^rQl{Q^^^-<Q^^^»-'Qf (46) 



(43) Using the F-matrix we calculate that 



(47) 

Prom figure 9 we then see that the effective operator be- 
tween a and b is 



Note that the two interactions must have the same /? by 
charge conservation (sec figure 10). The induced effective 
Hamiltonian between a and b at second order is 



(48) 



13 



where c is a constant containing the normaHzation of the 
projection operator T^^'') relative to its graphical repre- 
sentation, times the numerical factor one gets from re- 
ducing the portion of the graph in figure 10 between the 
two /3's using F-matrix moves. The product in front of 
comes out to 

Evaluating this expression using (34) we see that, up to 
/3-independent factors, it is equal to 

L2/3 + lJ-i/2 + L/3 + IJ^ - 1)-^ (50) 

Though this is not a decreasing function of /3, both of the 
factors are minimized in absolute value at (3=1, so the 



absolute value of the expression is maximized at /3 = 1. 
Thus we have shown that second order decimations also 
have the Q'-''-' as eigenvectors, and Q'^^^ has the highest 
magnitude eigenvalue. The second order decimation rule 
for /3 = 1 is then: 



H^^=^CjQ[% (51) 



where 



C, = \:^{l2j + l\,-l2j + l\-y\ (52) 
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This can be circumvented, however, by renormalizing the 
topological charges such that a certain fusion channel is 
eliminated. For instance, if we want to exclude the 1 ® 1 — »• 
2 channel between two sites, we can demote the two sites 
from charge 1 to 2, with the fusion channels now possible 
2 ® 2 = © 1 (this follows the procedure first suggested 
in*). In the limit of strong disorder, however, this is not 
necessary. 



